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I. INTRODUCTION 



There are several approaches to calculate the B — > light meson transition form factors, 
such as the lattice QCD technique, the QCD light-cone sum rules (LCSRs) and the per- 
turbative QCD (PQCD) approach. The PQCD calculation is more reliable when the in- 
volved energy scale is hard, i.e. in the large recoil regions; the lattice QCD results of the 
B — > light meson transition form factors are available only for soft regions; while, the QCD 
LCSRs can involve both the hard and the soft contributions below m\ — Im^x {x is a typical 
hadronic scale of roughly 500 MeV) and can be extrapolated to higher q 2 regions. Therefore, 
the results from the PQCD approach, the lattice QCD approach and the QCD LCSRs are 
complementary to each other, and by combining the results from these three methods, one 
may obtain a full understanding of the B — > light meson transition form factors in its whole 
physical region. In Refs.jl, 2], we have done a consistent analysis of the B — > ix transition 
form factor in the whole physical region. Similarly, one can obtain a deep understanding 
of the B K transition form factor in the physical energy regions by combining the QCD 
LCSR results with the PQCD results and by properly taking the SUf(3) breaking effects 
into account. 

The B — > K transition form factors are defined as follows: 

(K(p)\qj,b\B(p B )) = F^{q 2 ) ((p + M. ~ ^^gfr) +i^V) M| ~ 2 Ml % 
= 2F^ K (q 2 )p, + F B - K (q 2 )q,, (1) 

where the momentum transfer q = pb~P- If we confine ourselves to discuss the semi-leptonic 
decays B — > Klui, it is found that the form factors F B ~^ K (q 2 ) is irrelevant for light leptons 
(I = e, /x) and only F+~^ K (q 2 ) matters, i.e. 

^( B - A ^'-iiS A3/2(?2)|Frt ' tf)|2 - (2) 

where \{q 2 ) = (M§ + M\ — q 2 ) 2 — 4M^M^ is the usual phase-space factor. So in the 
following, we shall concentrate our attention on F+^ K (q 2 ). 

The B — > K transition form factor has been analyzed by several groups under the QCD 
LCSR approach , where some extra treatments to the correlation function either from 

the B-meson side or from the kaonic side are adopted to improve their LCSR estimations. It 
is found that the main uncertainties in estimation of the B —>■ K transition form factor come 
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from the different twist structures of the kaon wave functions. It has been found that by 
choosing proper chiral currents in the LCSR approach, the contributions from the pseudo- 
scalars' twist-3 structures to the form factor can be eliminated g, [7|. In the present paper, 
we calculate the B — > K form factor with chiral current in the LCSR approach to eliminate 
the most uncertain twist-3 light-cone functions' contributions. And more accurately, we 
calculate the 0(a s ) corrections to the kaonic twist- 2 terms. The S77/(3)-breaking effects 
from the twist-2 and twist-4 kaon wave functions shall also be discussed. 

In Reffl, we have calculated the B - K transition form factor up to 0(l/m?) in the 
large recoil region within the PQCD approach 8], where the B-meson wave functions \I/ b 
and ^> b that include the three- Fock states' contributions are adopted and the transverse 
momentum dependence for both the hard scattering part and the non-perturbative wave 
function, the Sudakov effects and the threshold effects are included to regulate the endpoint 
singularity and to derive a more reliable PQCD result. Further more, the contributions from 
different twist structures of the kaon wave function, including its S77/(3)-breaking effects, 
are discussed. So we shall adopt the PQCD results of Ref.Q] to do our discussion, i.e. to 
give a consistent analysis of the B — > K transition form factor in the large and intermediate 
energy regions with the help of the LCSR and the PQCD results. 

The paper is organized as follows. In Sec. II, we present the results for the B — > K 
transition form factor within the QCD LCSR approach. In Sec. Ill, we discuss the kaonic 
DAs with 577/(3) breaking effect being considered. Especially, we construct a model for the 
kaonic twist-2 wave function based on the two Gegenbauer moments af and ■ Numerical 
results is given in Sec. IV, where the uncertainties of the LCSR results and a consistent 
analysis of the B — > K transition form factor in the large and intermediate energy regions 
by combining the QCD LCSR result with the PQCD result is presented. The final section 
is reserved for a summary. 

II. Ff^ K (q 2 ) IN THE QCD LIGHT-CONE SUM RULE 

The sum rule for F?^ K (q 2 ) by including the perturbative 0(a s ) corrections to the kaonic 



m 

f B F^«(<j>) = £ e<»i->/" 2 [pg(s, ,f) + „£(o 2 )] ds , (3) 



twist-2 terms can be schematically written as 

1 r s o 
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where Pt2( s iQ 2 ) * s the contribution from the twist- 2 DA and Pj^iq 2 ) is f° r twist-4 DA, Jb 
is the B-meson decay constant. The Borel parameter M 2 and the continuum threshold sq 
are determined such that the resulting form factor does not depend too much on the precise 
values of these parameters; in addition the continuum contribution, that is the part of the 
dispersive integral from s to oo that has been subtracted from both sides of the equation, 
should not be too large, e.g. less than 30% of the total dispersive integral. The functions 
pif2 (s, q 2 ) and p^iq 2 ) can be obtained by calculating the following correlation function with 
chiral current 

n>,g) = i J rf 4 xe^(K(p)|T{s(x)7 At (l + 75)^)J(0)^(l + 75)rf(0)}|0> 

= n + [g 2 ,(p + g) 2 ]^ + n_[g 2 ,(p + g) 2 ] gAt . (4) 



The calculated procedure is the same as that of B — > tt form factor that has been done in 
Refs. , [10J ■ So for simplicity, we only list the main results for B — > K and highlight 

the parts that are different from the case of B — > tt, and the interesting reader may turn to 
Refs. 0,128 for more detailed calculation technology. 
As for p^ 2 {s, q 2 ), it can be further written as 



pT2(s,q 2 ) 



f. 



^ du(pK(u, /i)Im T T2 1 



TT JO 



ml 2 ' rat 2 ' 



u,n 



(5) 



where 7V 2 



1 s 

mZ* 1 m: 



b-quark pole mass 



W, u, p ) is the renormalized hard scattering amplitude, ml stands for the 



9j. Defining the dimensionless variables r\ = q 2 /ml 2 , r 2 = (p + q) 2 /m 
and p = [r\+ u(r 2 — ri) — u(l — ^M^/ml 2 }, up to order a s , we have 

lmT T2 {r 1 ,r 2 ,u,p,) 



*2 
b 



71 



*(i-,) + 2^foi-/>) 

Att 



7T 



6 + 3 In 



m 



*2 



P 1 



2Li 2 (n) 



2Li 2 ( 1 - / o) - 2 ( In ^ + 2 (\nr 2 + In ^ ' >~ 



\ — T\ 



r-2 



-9(p-l) 



Hp-K 



p-1 



1 1 



r 2 - p \p r 2 



+2 lnr 2 + 2 - 21n(r 2 - 1) + In 



m 



*2 * 



r 2 



p-1 



1 — Tx 



l-p 



2(1 -n) 



(ri - r 2 )(r 2 - p) 



ln^-21n P - 1 



r 2 



41np 

r 2 - 1 ) p-1 
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2(r 2 - i; 



(ri - r 2 )(p - n) 
9(1 -p) 



hip - 21n(p - 1) + 1 - In 



m 



*2 1 



+ 



2 lnr 2 + 21n(r 2 - 1) 

V r 2 



In 



*2 N 



2(1 -r 2 ) 2(1 -n) 

r 2 (r 2 -p) (ri - r 2 )(r 2 - p) 



1 + ln 



r 2 



p-1 



fa - 13 



In 



(6) 



for the case of r± < 1 and r 2 > 1. As for the coefficients of 5(1 — p), the higher power 
suppressed terms of order 0((M^/ml 2 ) 2 ) have been neglected due to its smallness. The 
dilogarithm function Li 2 (x) = — Jq y ln(l —t) and the operation " + " is defined by 



dpf(p): 



P 



dp[f{p) - /(!)] 



1-p 



(7) 



In the calculation, both the ultraviolet and the collinear divergences are regularized by 



dimensional regularization and are renormalized in the MS scheme with the totally anti- 
commuting 75. And similar to Ref. [3I, to calculate the renormalized hard scattering ampli- 



tude Tt2 



sy, it, pj, the current mass effects of s-quark are not considered due to their 
smallness. By setting Mx — > 0, it returns to the case of B — > ir and it can be found that 
the coefficients of 9(p — 1) and 9(1 — p) agree with those of Refsj7, 9], while the coefficients 
of 6(1 - p) confirm that ot Ref.fl and differ from that of Ref.ff The present results ean 
be checked with the help of the kernel of the Brodsky-Lepage evolution equation 11[ , since 
the p-dependences of the hard scattering amplitude and of the wave function should be 
compensate to each other. 

As for the sub- leading twist-4 contribution p^(q 2 ), we calculate it only in the zeroth 
order in a s , i.e. 



O b e^p%(q 2 )ds m ff K e 



Ml 



M 2 B 



due' 



f. (1 . u)(q 2_ uM 2 K) / 2 g 2 ( u ) 8m 2 [ gi (u) + G 2 (u)]" 



uM 2 



v?M A 



+ / dv Dai 



9(ai + va 3 



mj|-(l-a 1 -ua3)(q 2 -(a 1 +ua3)Mj^) 



x 



(ai + va 3 ) 2 M 2 
(2(p±(ai) + 2(p±(ai) - (p\\(ai) - <p[|(ati)jj 

where ip±(ai), ^j_(a»), <p\\(cxi) and {p\\(aii) are three-particle twist-4 DAs respectively, and 
gi(u) and g 2 (u) are two-particle twist-4 wave functions. Here, G 2 (it) = Jq g 2 (v)dv, A = 



(8) 



v /(.o-g 2 -A4)2+4Afg.(m2-g2)_ (so _ g 2_ M | ) 



2M% 



and Sq denotes the subtraction of the continuum from 



the spectral integral. By setting M% — > (the lower integration range of u should be changed 



to be A 



S0~<2 



for the case), we return to the results of B — > n 
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III. THE DISTRIBUTION AMPLITUDES OF KAON 



twist- 2 DA moments 



Generally, the leading twist-2 DA (ftx can be expanded as Gegenbauer polynomials: 



4> K (u,n ) = 6u(l - u) 



l + E«n(/^ /2 (2« 



n=l 



(9) 



In the literature, only af(fi Q ) is determined with more confidence level and the higher 
Gegenbauer moments are still with large uncertainty and are determined with large errors. 
Alterative determinations of Gegenbauer moments rely on the analysis of experimental data. 
The first Gegenbauer moment a? has been studied by the light-front quar k model 12), 

the LCSR approach B Q, H, y , Q 



Rcf. 



and the lattice calculation 



18 



19| and etc. In 



14j , the QCD sum rule for the diagonal correlation function of local and nonlocal 



axial-vector currents is used, in which the contributions of condensates up to dimension 
six and the (9(a s )-corrections to the quark-condensate term are taken into account. The 
moments derived there are close to that of the lattice calculation 18|, [l9| , so we shall take 
af(lGeV) = 0.05 ± 0.02 to do our discussion. At the scale p b = J M% - mf ~ 2.2 GeV, 



af (> 6 ) = 0.793af (lGeV) with the help of the QCD evolution. 



The higher Gegenbauer moments, such as a^, are still determined with large uncertainty 



and are determined with large errors 



For example, Ref. [2lj shows that 



the value of aff is very close to the asymptotic distribution amplitude, i.e. \a^(lGeV )\ < 



0.04; while Refs. |20j ] gives larger values for af , i.e. af ( lGeV) = 0.16 ± 0.10 



(lGeV) = 0.27lg; 



and af(2GeV) = 0.175 ± 0.065 |20j. It should be noted that 



the value of aff affects not only the twist-2 structure's contribution but also the twist-4 
structures' contributions, since the 5't/j-(3)-breaking twist-4 DAs also depend on af" due to 
the correlations among the twist-2 and twist-4 DAs as will be shown in the next subsection. 
Since the value of af" can not be definitely known, we take its center value to be a smaller 
one, i.e. aff^GeV) = 0.115, for easily comparing with the results of Ref. [3|. Further more, 
to study the uncertainties caused by the second Gegenbauer moment aff , we shall vary a 



K 



within a broader region, e.g. a? (lGeV) G [0.05,0.15], so as to see which value is more 



favorable for a* by comparing with the PQCD results. 
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B. Models for the twist-2 and twist-4 DAs 



Before doing the numerical calculation, we need to know the detail forms for the kaon 
twist-2 DA and the twist-4 DAs. 

As for the twist-2 DA, we do not adopt the Gegenbauer expansion Q, since its higher 
Gegenbauer moments are still determined with large errors whose contributions may not be 
too small, i.e. their contributions are comparable to that of the higher twist structures. For 
example, by taking a typical value af(lGeV) = —0.015 our numerical calculation shows 
that its absolute contributions to the form factor is around 1% in the whole allowable energy 
region, which is comparable to the twist-4 structures' contributions. Recently, a reasonable 
phenomenological model for the kaon wave function has been suggested in Ref.fl, which is 
determined by its first Gegenbauer moment af, by the constraint over the average value of 



the transverse momentum square, (k^)]^ 2 0.350GeV [22j, and by its overall normalization 
condition. With the help of such model, a more reliable PQCD calculation on the B —>■ K 
transition form factors up to 0(l/ml) have been finished. 

In the following, we construct a kaon twist-2 wave function following the same arguments 
as that of Ref.Q] but with slight change to include the second Gegenbauer moment 's 
effect, i.e. 



V K (x,k ± ) = \l+B K Cl ,2 {2x-l)+C K C ? f{2x-l)} .f K . exp 

x 1 — X) 



X 1 — X ) 

(10) 

where q = u, d, Cf^l — 2x) are the Gegenbauer polynomial. The constitute quark masses 
are set to be: m q = 0.30GeV and m s = 0.45GeV. The four parameters A K , B K , Ck and 



Pk can be determined 



yy the first two Gegenbauer moments af and , the constraint 



(kj_)i/ 2 « 0.350GeV 22] and the normalization condition Jq dx $k\<nl j^-^k{x, kj_) = 1. 
For example, we have A K (/j, b ) = 252.044GeK" 2 , B K (fx b ) = 0.09205, C K (fi b ) = 0.05250 and 
p K = 0.8657Ge^- 1 for the case of af(lGeV) = 0.05 and af {iGeV) = 0.115. Quantita- 
tively, it can be found that Bk, Ck and (3k decreases with the increment of af^; (3k decreases 
with the increment of , while Bk and Ck increase with the increment of ■ Under such 
model, the uncertainty of the twist-2 DA mainly comes from af" and ■ It can be found that 
the SUf(3) symmetry is broken by a non-zero Bk and by the mass difference between the s 
quark and u (or d) quark in the exponential factor. The SUf(3) symmetry breaking effect 



of the leading twist kaon distribution amplitude has been studied in Refs. 



14 



23j and refer- 



ences therein. The SUf(3) symmetry breaking in the lepton decays of heavy pseudoscalar 



mesons and in the semileptonic decays of mesons have been studied in Ref. 



2J] . After doing 



the integration over the transverse momentum dependence, we obtain the twist- 2 kaon DA, 



4> K (x,/i ) 



d 2 k. 



k 2 ± <fl% 167T 3 
A K 



16vr 2 /3 2 
x exp 



^ x (x,kj 



1 + B K C 3 1 /2 (2x - 1) + C K Ct'\2x - 1 



3/2, 



-Pk — + 



mZ 



x 



1 — x 



exp 



x(l — x) 



(11) 



where /xq = A*6 for the present case. Then, the Gegenbauer moments a^(/io) can be defined 



as 



a n (/i ) = 

where </>#(! — x, /i ) other than 



/ 1 cia: ( /»x(l-x,/io)C^ 2 (2x-i; 
^^(l -x)[C^ /2 (2x- l)] 2 



(12) 



defined in the literature, e.g. 



13 



14 



x, lip) is adopted to compare the moments with those 



15j, since in these references x stands for the momentum 



fraction of s-quark in the kaon (K), while in the present paper we take x as the momentum 
fraction of the light g-(anti) quark in the kaon (K). 

The twist-3 contribution is eliminated by taking proper chiral currents under the LCSR 
approach, so we only need to calculate the subleading twist-4 contributions. The needed 



four three-particle twist-4 DAs that are defined in Ref. 



251 ] can be expressed as |2i 



(p±(oa) = 30a 2 (a 2 -a:i' 



h 00 + V«3 + ^ ^io(5« 3 - 3) 



<P\\(<*i) 

where 



-30e^ 



h 00 (l - a 3 ) + hoi [«3(1 - "3) 
120aia 2 a 3 [a 10 (ai - a 2 )} , 
120aia 2 tt3 boo + ^io(3a 3 - 1)] , 



6ai«2 



+ h 



10 



a 3 (l - a 3 ) - -(a 2 + a 2 2 ) 



(13) 



^-00 


Ml 
= v 00 = 3 


7/4 = " 


a w 


21M 2 
g V^i ■ 


9 

"20 a 


v w 


21M 2 K 




h i 


7M 2 K 
= 4 


20 2 



' 3 



2MI 



9 

20 



S 2 e - ^agMjc, 



K 



l S 2 e - -a«M 2 
3 20 2 



K 



1 Similar to Ref.Q, we adopt the results that only include the dominant meson-mass corrections. The less 
important meson-mass correction terms are not taken into consideration. 
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and 

ho = + ±a*M% = -/e + laf M* K , 

with rj 4 = 5 2 /M 2 K , u 4 = 8e/21 and 5 2 (lGeV) = 0.20GeV 2 and e(lCfey) = 0.53 |26j. With 
the help of QCD evolution, we obtain <5 2 (u{,) = 0.16GeV 2 and £(//&) = 0.34. It can be found 
that the dominant meson-mass effect are proportional to a 2 and MjL so if setting — > 



or the value of is quite small, then we return to the results of Ref . 25j] . For the remaining 
two-particle twist-4 wave functions, their contributions are quite small in comparison to the 
leading twist contribution and even to compare with those of the three-particle twist-4 wave 
functions. And by taking the leading meson-mass effect into consideration only, they can be 
related to the three-particle twist-4 wave functions through the following way: 

ru ru j 

g 2 (u) = doix I da 2 — [2ip±(ai) - wdan)] (14) 
Jo Jo a 3 



and 



g 1 (u)+ I dvg 2 (v) = ^ / da x j da 2 -^(ua 1 - ua 2 )[2ip ± (ai) - (p\\{(Xi)], (15) 



o ' 2 Jo Jo a 



which lead to 



mm 

gi(u) = — [— 5uu(9h 00 + 3h 01 — 6hi + 4uho\U + Wuhi u) + a w (6 + uu(9 + 80mm))] + 
6 

ai M 3 (10 - 15m + 6m 2 ) In it + ai M 3 (10 - 15m + 6m 2 ) \nu, (16) 
huuiu — u) 

g 2 {u) = [Ah 00 + 8ai uu - h 10 (l + 5mm) + 2/i i(l - mm)] . (17) 

Similarly, it can be found that when setting af — > 0, the above expressions of g\{u) 



and g 2 {u) return to those of Ref.[25j. Here by adopting the relations 4- g 2 



u 



-|lim M 2,^ M^[g K (u) - <j) K (u)] and gi(u) - Jq dvg 2 (v) = ± lini M a_> M%A(u), one can 
conveniently obtain the higher mass-correction terms for g\{u) and g 2 (u) on the basis of 



gx{u) and A(m) derived in Refs. 13|, l26j, and numerically, it can be found that these terms' 
contributions are indeed small. 

IV. NUMERICAL RESULTS 
A. basic input 

In the numerical calculations, we use 

M B = 5.279GeV, M K = A9AMeV, f K = 160MeV, / w = 131MeV. (18) 



Total 

Twist-2 

Twist-4 



M 2 (GeV 2 ) 



FIG. 1: Ff^ K (q 2 ) as a function of Borel parameter M 2 at q 2 = 6GeV 2 , where s = 33.5GeV 2 , 
af (lGeV) = 0.05, af (lGeV) = 0.115, ml = A.7GeV. The solid line stands for the total contribu- 
tions, the dashed line is for NLO result of the twist-2 kaonic wave function and the dotted line is 
for the LO result of twist-4 kaonic wave functions. 

Next, let us choose the input parameters entering into the QCD sum rule. In general, 
the value of the continuum threshold sq might be different from the phenomenological value 
of the first radial excitation mass. Here we set the threshold value of sq to be smaller 
than s™ ax ~ 34GeV 2 , whose root is slightly bigger than the mass of the B-meson first 
radlal I xcitation predicted by t* potent m ode> 0. The po,e qUark mass ml is taken 
as 4.7 — 4.9GeV. Another important input is the decay constant of B meson To keep 
consistently with the next-to-leading order calculation of twist-2 contribution, we need to 
calculate the two-point sum rule for fg up to the corrections of order a s . And in doing the 
numerical calculation, we shall adopt the NLO fs to calculate the NLO twist-2 contribution 
and LO fs for the LO twist-4 contributions for consistence. 

The reasonable range for the Borel parameter M 2 is determined by the requirement that 
the contributions of twist-4 wave functions do not exceed 10% and those of the continuum 
states are not too large, i.e. less than 30% of the total dispersive integration. At a typical 
q 2 = 6GeV 2 , we draw F+^ K (q 2 ) versus M 2 in Fig.([T]). It can be found that the contribution 
from the kaonic twist-2 wave function slightly increases with the increment of M 2 while the 
contributions from the kaonic twist-4 wave functions decreases with the increment of M 2 , 
as a result, there is a platform for F^ K (q 2 ) as a function of the Borel parameter M 2 for 
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the range of 8GeV 2 < M 2 < 18GeV 2 . For convenience, we shall always take M 2 = 12GeV 2 
to do our following discussions. 



B. uncertainties for the LCSR results 

In the following we discuss the main uncertainties caused by the present LCSR approach 
with the chiral current. 

The present adopted chiral current approach has a striking advantage that the twist-3 
light-cone functions which are not known as well as the twist-2 light-cone functions are 
eliminated, and then it is supposed to provide results with less uncertainties. In fact, it 
has been pointed out that the twist-3 contributions can contribute ~ 30 — 40% to the total 



contribution [28J by using the standard weak current in the correlator, e.g. 

IL>,g) = i J d i xe^%K{p)\T{s{x) 1 ^b{x)X^d^W > • (19) 

If the twist-3 wave functions are not known well, then the uncertainties shall be large 2 . 
So in the literature, two ways are adopted to improve the QCD sum rule estimation on 
the twist-3 contribution: one is to calculate the above correlator by including one-loop 
radiative corrections to the twist-3 contribution together with the updated twist-3 wave 
functions [3!]; the other is to introduce proper chiral current into the correlator, cf. Eq.(j4]), 
so as to eliminate the twist-3 contribution exactly, which is what we have adopted. We 



shall make a comparison of these two approaches in t 
adopt the following form for the QCD sum rule of Ref. 
contributions from different Gegenbauer moments: 



re following. For such purpose, we 
3J], which splits the form factor into 



F^ K (q 2 ) = f as (q 2 ) + af (fi )f< (q 2 ) + af (n )F> (q 2 ) + af(n )f< (q 2 ), (20) 

where f as contains the contributions to the form factor from the asymptotic DA and all 
higher- twist effects from three-particle quark-quark-gluon matrix elements, / 0l '° 2 '° 4 con- 
tains the contribution from the higher Gegenbauer term of DA that is proportional to af, 
af and af respectively. The explicit expressions of f as > a i > a 2 > a A can be found in Table V 



2 A better behaved twist-3 wave function is helpful to improve the estimations, e.g. Ref. [29j provides such 
an example for the pionic case. 
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0.0 ■ 



I 

10 



15 



q 2 (GeV 2 ) 



FIG. 2: Ff^ K {q 2 ) for af{lGeV) G [0.03,0.07], a$(lGeV) G [0.05,0.15] and m* G [4.7,4.9]GeV. 
The solid line is obtained with of (lGeV) = 0.03, of (lGeV) = 0.15 and m£ = 4.9GeF; the dashed 
line is obtained with of(lGeV) = 0.07, af(lGeF) = 0.05 and raj* = 4.7GeF, which set the upper 
and the lower ranges of F^?~* (q 2 ) respectively. As a comparison, the shaded band shows the result 
of Ref.[3] together with its 12% theoretical uncertainty. 



and Table IX of Ref. [3|. And in doing the comparison, we shall take the same DA moments 
for both methods, especially the value of af (//o) is determined from Eq. (TT2]) . 

We show a comparison of our result of F^~* K (q 2 ) with that of Eq. fl20l in Fig. (j2J) by varying 
af(lGeV) G [0.03,0.07], af(lGeV) e [0.05,0.15] and m* b G [4.7,4.9]GeV. In Fig.® the 
solid line is obtained with af(lGeV) = 0.03, af (lGeV) = 0.15 and m* = 4.9Ge\/; the 
dashed line is obtained with af(lGeV) = 0.07, af (lGeV) = 0.05 and m£ = 4.7GW, which 
set the upper and the lower ranges of F+^ K (q 2 ) respectively. The shaded band in the figure 
shows the result of Eq. fl20|) within the same af and af region and with its 12% theoretical 
uncertainty [3j|. It can be found that our present LCSR results are consistent with those of 
Ref. (3]] within large energy region q 2 G [0, IhGeV 2 ]. In another words these two treatments 
on the most uncertain twist-3 contributions are equivalent to each other, while the chiral 
current approach is simpler due to the elimination of the twist-3 contributions. One may 
also observe that in the lower q 2 region, different from Ref.[3j where F^ K (q 2 ) increases with 
the increment of both af and af , the predicted F+^ K (q 2 ) will increase with the increment 
of af but with the decrement of af . This difference is caused by the fact that we adopt 
the model wave function fllOp to do our discussion, whose parameters are determined by the 



12 





LO result 


NLO result 


_ 


so 


M 2 


fB 


so 


M 2 


fB 


nib = 4.7 


33.5 


2.80 


0.165 


33.5 


2.80 


0.219 


m h = 4.8 


33.2 


2.39 


0.131 


33.2 


2.31 


0.174 


m h = 4.9 


32.8 


2.16 


0.0997 


32.8 


2.02 


0.132 



TABLE I: Parameters for where nib and /b are given in GeV , so and M in GeV . 

combined effects of af and af^ while in Ref.js], af" and are varied independently and 
then their contributions are changed separately. 

Next we discuss the main uncertainties caused by the present LCSR approach with the 
chiral current. Firstly, we discuss the uncertainties of F^ K (q 2 ) caused by the effective 
quark mass ml by fixing af{lGeV) = 0.05GeV and (lGeV) = O.llSGeV^. Under such 
case, the value of sq, the LO and NLO vales of fs should be varied accordingly and be 
determined by using the two-point sum rule with the chiral currents, e.g. to calculate the 
following two-point correlator: 

Tl(q 2 )=iJ d 4 xe^(0|g(x)(l + 75 )^)J(0)(l-75)g(0)|0). (21) 

The sum rule for fs up to NLO can be obtained from Ref. through a proper combination 
of the scalar and pseudo-scalar results shown there 3 , which can be schematically written as 



f B M B e 



-Ml/M 2 



p tot {s)e- s ' h 1 ds, 



(22) 



where the spectral density p tot (s) can be read from Ref. 30|. The Borel parameter M 2 and 
the continuum threshold So are determined such that the resulting form factor does not 
depend too much on the precise values of these parameters; in addition, 1) the continuum 
contribution, that is the part of the dispersive integral from so to oo, should not be too large, 
e.g. less than 30% of the total dispersive integral; 2) the contributions from the dimension- 
six condensate terms shall not exceed 15% for fg. Further more, we adopt an extra criteria 
as suggested in Ref. [3[] to derive fs- i.e. the derivative of the logarithm of Eq. (l22l with 



3 One needs to change the c-quark mass to the present case of 6-quark mass and we take (^f G^„G aAlI/ ) 
2 x (0.33Gey) 4 [EiJ and a s (qq} 2 = 0.162 x lCT 3 Gey 6 to do the numerical calculation. 
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FIG. 3: F+^ K (q 2 ) as a function of q 2 with varying m£. The solid, dashed line and the dash-dot 
line are for m* b = 4.7GeV, 4.8GeV and 4.9GeV respectively, where af{lGeV) = 0.05GeF and 
af(lGeV) = 0.115GeV. 

respect to 1/M 2 gives the B-meson mass Mb, 

Ml = r p tot (s)e- s ' M2 sds I f° p tot (s)e- s ' M2 ds, 



b= ,p y s > e sas / / , p 

and we require its value to be full-filled with high accuracy ~ 0.1%. These criteria define 
a set of parameters for each value of m£. Some typical values of fs are shown in TAB JH 
where fs is taken as the extremum within the reasonable region of (M 2 ,sq) and the value 



of ml is taken as 32|: ml ~ 4.8 ± O.lGeV . fs decreases with the increment of m£. The 
NLO result agrees with the first direct measurement of this quantity by Belle experiment 
Jb = 229^3 < j i (stat)^3y(syst) MeV from the measurement of the decay B~ — > rv T {33 ]. 

The value of F?~* K (q 2 ) for three typical values of mj, i.e. ml = 4.7GeV, 4.8GeV and 
4.9GeV respectively, are shown in Fig.Q. F^~* K (q 2 ) increases with the increment of ra£. It 
can be found that the uncertainty of the form factor caused by ml G [0A7GeV, 0A9GeV] is 



5% at q 2 = and increases to ~ 9% at q 2 = 14GeV 2 . By taking a more accurate 



m b , e.{ 



ml = (4.80 ± 0.0B)GeV as suggested by Ref. [3j, the uncertainties can be reduced to ~ 3% 
at q 2 = and ~ 5% at q 2 = UGeV 2 . 

Secondly, we discuss the uncertainties of F^ K (q 2 ) caused by the twist-2 wave function 
^k, i-e. the two Gegenbauer moments af (lGeV) and af (lGeV). For such purpose, we 
fix s = 33.5GeV 2 and ml = 4.7GeV. To discuss the uncertainties caused by af (lGeV), 
we take a%(lGeV) = 0.115. F^ K (q 2 ) for three typical af (lGeV), i.e. af (lGeV) = 0.03, 
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FIG. 4: F^ K {q 2 ) as a function of q 2 with varying of (lGeV), where a£(lGeV) = 0.115. The solid 
line, the dashed line and the dash-dot line are for of (lGeV) = 0.03, 0.05 and 0.07 respectively. 




FIG. 5: F^ K (q 2 ) as a function of q 2 with varying of (lGeV), where of (lGeV) = 0.05. The solid 
line, the dashed line and the dash-dot line are for of (lGeV) = 0.05, 0.10 and 0.15 respectively. 

0.05 and 0.07 respectively, are shown in Fig.(j4]). F+^ K (q 2 ) decreases with the increment of 
af . It can be found that the uncertainty of form factor caused by af (lGeV) G [0.03,0.07] 
is small, i.e. it is about 3% at q 2 = and becomes even smaller for larger q 2 . Similarly, 
to discuss the uncertainties caused by af- (lGeV), we fix af (lGeV) = 0.05. Since the 
value of is less certain than af , so we take three typical values of af (lGeV) with 
broader separation to calculate F^ K (q 2 ), i.e. af (lGeV) = 0.05, 0.10 and 0.15 respectively. 
The results are shown in Fig. ([5]). It can be found that the uncertainty of the form factor 
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caused by af (lGeV) G [0.05,0.15] is also small, i.e. it is about 5% at q 2 = and becomes 
smaller for larger q 2 . F+^ K (q 2 ) increases with the increment of af in the lower energy 
region q 2 < lOGeV 2 and decreases with the increment of af in the higher energy region 
q 2 > lOGeV 2 . 

As a summary, a more accurate values for ml, af and af shall be helpful to derive a more 
accurate result for the form factor. Our results favor a smaller af to compare with the form 
factor in the literature, e.g. af (lGeV) < 0.15. And under such region, the uncertainties 
from af is small, i.e. its uncertainty is less than 5% for af (lGeV) G [0.05,0.15]. It can 
be found that by varying af (lGeV) G [0.03,0.07] and af(lGeV) G [0.05,0.15], the kaonic 
twist-4 wave functions' contribution is about 6% of the total contribution at q 2 = 0. The 
uncertainties of af shows that the S'L r j(3)-breaking effect is small but it is comparable 
to that of the higher twist structures' contribution. So the SUf(3) breaking effect and 
the higher twist's contributions should be treated on the equal footing. Using the chiral 
current in the correlator, as shown in Eq.(j4|), the theoretical uncertainty can be remarkably 
reduced. And our present LCSR results are consistent with those of Ref. 3]] within large 
energy region q 2 G [0, 15GeV 2 ], which is calculated with the correlator ffTUj) and includes 
one-loop radiative corrections to twist-2 and twist-3 contributions together with the updated 
twist-3 wave functions. In another words these two approaches are equivalent to each other 
in some sense, while the chiral current approach is simpler due to the elimination of the 
more or less uncertain twist-3 contributions. For higher energy region q 2 > 15GeV 2 , the 
LCSR approach is no longer reliable. Therefore the lattice calculations, would be extremely 
useful to derive a more reliable estimation on the high energy behaviors of the form factors. 



C. SUf(3) breaking effect of the form factor within the LCSR 

To have an overall estimation of the SUf(3) breaking effect, we make a comparison 
of the B —* 7T and B — > K form factors: F^ n (q 2 ) and F+^ K \q 2 ). The formulae 
for F+^ n (q 2 ) can be conveniently obtained from that of F+^ K (q 2 ) by taking the limit 
Mk — > 0. In doing the calculation for F+~* n (q 2 ), we directly use the Gegenbauer ex- 
pansion for pion twist-2 DA, because different to the kaonic case, now the higher Gegen- 
bauer terms' contributions are quite small even in comparison to the twist-4 contributions, 
e.g. by taking aJ(lGeV) = —0.015 [3], our numerical calculation shows that its abso- 
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FIG. 6: Comparison of F^ K \q 2 ) and Ff ^(q 2 ), where m* b = 4.7GeV, s = 33.5GeV 2 , /|° = 
0.l65GeV, /^ LO = 0.219GeV, of (lGeV) = 0.05 and a 2 /K (lGeV) = 0.115. The solid line and 
the dashed line are for F?^ K (q 2 ) and FP~* n (q 2 ) respectively. 



lute contributions to the form factor is less than 0.5% in the whole allowable energy re- 
gion. We show a comparison of F^ K (q 2 ) and F^ 1T (q 2 ) in Fig. with the parame- 
ters taken to be m* b = A.7GeV, s = 33.5GeV 2 } f%° = 0.165GeV, fg LO = 0.219GeV, 
af(lGeV) = 0.05 and a 2 /K (lGeV) = 0.115. Secondly, by varying m* b E [4.7,4.9]GeV, 
a?(lGeV) E [0.03,0.07] and a 2 /K (\GeV) E [0.05,0.15], we obtain Ff^(0) E [0.239,0.294] 
and Ff^ A '(0) E [0.273,0.349]. Then we obtain ^g^M = 1.16 ± 0.03, which favors a small 



r (o) 



SUf(3) breaking effect and is consistent with the PQCD estimation 1.13 ± 0.02 
QCD sum rule estimations, e.g. [Ff -* K (0)/F^(0)] « 1.16 [3] 4 , 1.08±g; 19 



.17 



the 



34| and 1.36+ 012 



-0.09 



14| respectively, and a recently relativistic treatment that is based on the study of the 



Dyson-Schwinger equations in QCD, i.e. [Ff ^(0)/Ff ^(0)] = 1.23 35]. 



D. consistent analysis of the form factor within the large and the intermediate 
energy regions 

Recently, Ref 8j gives a calculation of the B — > K transition form factor up to 0(1/ m 2 ) 
in the large recoil region within the PQCD approach 8], where the B-meson wave functions 



To estimate the ratio [F?^ K (0)/Ff ^(0)] from Ref. 13], we take af(\GeV) = 0.05 ± 0.02 
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FIG. 7: LCSR and PQCD results for F^ K (q 2 ). The solid line is for our LCSR result, the dash-dot 
line is for the LCSR result of Ref.[3] with of (lGeV) = 0.07 and ag(lGeV) = 0.05. The shaded 
band is the PQCD result with A £ [0.50, 0.55] and 5 G [0.25, 0.30], where the dashed line is for the 
center values A = 0.525 and 5 = 0.275, the upper edge of the band is for A = 0.50 and 5 = 0.30 
and the lower edge of the band is for A = 0.55 and 5 = 0.25. 

\& b and b that include the three- Fock states' contributions are adopted and the transverse 
momentum dependence for both the hard scattering part and the non-perturbative wave 
function, the Sudakov effects and the threshold effects are included to regulate the endpoint 
singularity and to derive a more reliable PQCD result. Further more, the uncertainties for 
the PQCD calculation of the B —> K transition form factor has been carefully studied in 
Ref.[8|]. So we shall adopt the PQCD results of Ref. [8J] to do our discussion. Only we need 
to change the twist-2 kaon wave function \I/ k used there to the present one as shown in 
Eq.flTDD. 

We show the LCSR results together with the PQCD results in Fig. (J7J) . In drawing 
the figure, we take af (lGeV) = 0.07, af {IGeV) = 0.05 and m* b = A.8GeV. And the 
uncertainties of these parameters cause about ~ 10% errors for the LCSR calculation. While 
for the PQCD results, we should also consider the uncertainties from the B-meson wave 
functions, i.e. the values of the two typical parameters A and 5, and we take A G [0.50, 0.55] 
and S G [0.25,0.30] 8]. It can be found that the PQCD results can match with the LCSR 
results for small q 2 region, e.g. q 2 < lOGeV 2 . Then by combining the PQCD results with 
the LCSR results, we can obtain a consistent analysis of the form factor within the large and 
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the intermediate energy regions. Inversely, if the PQCD approach must be consistent with 
the LCSR approach, then we can obtain some constraints to the undetermined parameters 
within both approaches. For example, according to the QCD LCSR calculation, the form 
factor F^ K (q 2 ) increases with the increment of b-quark mass, then the value of m& can not 
be too large or too small 5 , i.e. if allowing the discrepancy between the LCSR result and the 
PQCD results to be less than 15%, then m£ should be around the value of 4.8 ± O.lGeV. 

V. SUMMARY 

In the paper, we have calculated the B —>■ K transition form factor by using the chiral 
current approach under the LCSR framework, where the SUf(3) breaking effects have been 
considered and the twist-2 contribution is calculated up to next-to-leading order. It is found 
that our present LCSR results are consistent with those of Ref. |3[ within large energy region 
q 2 G [0, 15GeV 2 }, which is calculated with the conventional correlator ( Tl9|) and includes one- 
loop radiative corrections to twist-2 and twist-3 contributions together with the updated 
twist-3 wave functions. And our present adopted LCSR approach with the chiral current is 
simpler due to the elimination of the more or less uncertain twist-3 contributions. 

The uncertainties of the LCSR approach have been discussed, especially we have found 
that the second Gegenbauer moment af prefers asymptotic-like smaller values. By varying 
the parameters within the reasonable regions: m% G [4.7, 4.9]GeV, af (lGeV) G [0.03,0.07] 
and a 2 /K (\GeV) G [0.05,0.15], we obtain Ff^(0) = 0.267± 0.026 and F^ K (0) = 0.311 ± 
0.038, which are consistent with the PQCD and the QCD sum rule estimations in the 

F B ^ K (0) 

literature. Consequently, we obtain pu^n^ = 1-16 ± 0.03, which favors a small SUf(3) 
breaking effect. Also, it has been shown that one can do a consistent analysis of the B — > K 
transition form factor in the large and intermediate energy regions by combining the QCD 
LCSR result with the PQCD result. The PQCD approach can be applied to calculate the 
B —>■ K transition form factor in the large recoil regions; while the QCD LCSR can be 
applied to intermediate energy regions. Combining the PQCD results with the QCD LCSR, 
we can give a reasonable explanation for the form factor in the low and intermediate energy 
regions. Further more, the lattice estimation shall help to understand the form factors' 

5 Another restriction on mt, is from the experimental value [i^ l on fg. 
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behaviors in even higher momentum transfer regions, e.g. q 2 > 15GeV 2 . So, we suggest 
such a lattice calculation can be helpful. Then by comparing the results of these three 
approaches, the B — > K transition form factor can be determined in the whole kinematic 
regions. 
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